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Abstract 

We present a fully non-perturbative computation of the mass of the b-quark in the 
quenched approximation. Our strategy starts from the matching of HQET to QCD in 
a finite volume and finally relates the quark mass to the spin averaged mass of the B s 
meson in HQET. All steps include the terms of order A 2 /mb- Expanding on [1], we 
discuss the computation and renormalization of correlation functions at order l/m^. 
With the strange quark mass fixed from the Kaon mass and the QCD scale set through 
tq = 0.5 fm, we obtain a renormalization group invariant mass M\, = 6.758(86) GeV or 
mb(fflb) = 4.347(48)GeV in the MS scheme. The uncertainty in the computed A 2 /mb 
terms contributes little to the total error and A 3 /m^ terms are negligible. The strategy 
is promising for full QCD as well as for other B-physics observables. 
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1 Introduction 



The mass of the b-quark, Mb, is a relevant input parameter for phenomenology analysis 
based on perturbation theory. Let us just mention the extraction of V n \> from inclusive 
b-decays [2, 3]. Mb is a fundamental parameter of the Standard Model of particle 
physics. Thus it should be determined precisely. One may of course turn the very first 
observation around. For instance, applying high order perturbation theory to sufficiently 
well integrated cross sections, the quark mass can be determined [4-15]. Still, the 
achievable precision is limited by the intrinsic uncertainty of perturbation theory and 
maybe more by experimental difficulties. On the other hand, the use of lattice QCD 
offers a strategy to compute the fundamental renormalization group invariant (RGI) 
parameters of QCD with very precisely known experimental input, e.g. N{ (= number 
of quarks flavours) meson masses as well as the nucleon mass; see [16] for a basic 
introduction. However, for the b-quark mass, such a computation is more involved than 
for the light quarks because the achievable inverse lattice spacings are below the mass of 
the quark. Effective theories have to be employed in a numerical treatment of the bound 
states. The most serious problem that arises is that a power law divergent (~ go/a) 
additive renormalization of the mass is present due to the absence of a chiral symmetry 
in the effective theories (even in the continuum) . Although at the lowest order in Heavy 
Quark Effective Theory (HQET) the subtraction is known to order g^/a [17-19], an (in 
the continuum limit) divergent remainder is unavoidable and the total uncertainty is 
difficult to estimate as long as the renormalization is carried out perturbatively. 

In [1] a general strategy was described which allows HQET at zero velocity to be 
implemented non-perturbatively on the lattice, including all renormalizations. 



experiment 



Lattice with am. 



m B = 5.4GeV $i(Li, M b ), $ 2 (ii, M b ) 

L2 = 2Li 



$? QET (L 2 ), $^ ET (i 2 ) , gm(Ml) *? QET (Li), ^^(Li) 

af"( Wl ),4 in ( Ul ) 

The basic idea is illustrated in the above diagram. It is founded on the knowledge 
of the relation between the RGI mass and the bare mass in QCD [20,21]. In a finite 
volume of extent L\ ~ 0.4 fm, one chooses lattice spacings a sufficiently smaller than 
1 / mb , such that the b-quark propagates correctly up to controllable discretization errors 
of order a 2 . Finite volume observables <E>j(Li,Mb) may then be computed as a function 
of the RGI mass Mb including an extrapolation to the continuum limit. The resulting 
values are equated to their representation in HQET - a step called matching, indicated 
by the r.h.s. of the diagram. Choosing now L\/a = O(10), with the same physical value 
of Li, one uses the knowledge of &i(L\, M^) to determine the bare parameters in the 
effective theory for a-values of about 0.025 fm to 0.05 fm. At these lattice spacings one 
then computes the same observables in a larger volume L2 = 2L\. Again these observ- 
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ables can be extrapolated to the continuum limit. Next, the knowledge of $j(L2,Mb) 
and the choice Lija = O(10) yields the bare parameters of the effective theory for a 
around 0.05 fm to 0.1 fm. One then has full control over the effective theory at lattice 
spacings where large volume observables, such as the B-meson mass, can be computed. 
Perturbation theory is completely avoided with the power divergent subtractions being 
taken care of non-perturbatively. 

We return to the specific application of computing Mb . The whole chain allows to 
express m-Q in terms of <J>j(Li, Mb) and thus as a function of Mb. This function naturally 
splits into various pieces which may be computed individually as they separately have a 
continuum limit. In particular, the step scaling functions a relate <f>j(Li) to §i{L2). As 
we will see below, at first order in 1/mb, two matching observables $X; ^2 are sufficient 
if we consider the spin averaged B-meson mass. 

The strategy requires all considered observables to be accurately described by the 

1 / ?7ib expansion. Naive counting estimates the accuracy of the quark mass as A x A^- 

and A x 2 2 ■ For a typical QCD scale A ~ 400 MeV both these terms yield the 

L 1 m b 

same, very small, estimate. In [22] the l/m-b expansion was tested for an even smaller 
L = Lq = Li/2 and found to be well behaved, as it is also the case in perturbation 
theory [23]. Here we will have additional cross checks by choosing different quantities 
$j in the matching step. 

In Sect. [2] we will go through the definition of the effective theory in order to fix 
some notations and give rules how the 1/mb expansion is implemented in practice. We 
also discuss correlation functions in the Schrodinger functional [24,25], which defines 
our finite volume geometry. These correlation functions are then used in Sect. [3] to form 
suitable dimensionless observables followed by a section which lists the step scaling 
functions. Sect. U] discusses the final formula for the RGI b-quark mass Mb. Numerical 
results for all quantities in the quenched approximation are discussed in Sect. [5j This 
includes also results from an alternative strategy as a check on the smallness of the 
1/m^ terms. 

2 Heavy quark effective theory on the lattice 

We start from the Eichten Hill static quark Lagrangian [26], using the notation of 
[27], but setting the mass counter term 5m to zero. Its effect is taken into account 
in the overall energy shift m-bare between the effective theory and QCD. Thus rribare is 
regularization dependent with a ~ 9q/cl divergence. For the sake of a light notation, 
we also drop the superscript W [27] for the different lattice discretizations of the static 
Lagrangian, but in the numerical computations these different versions will be used and 
referred to exactly as in that reference. We remind the reader that they differ only by 
the choice of the covariant derivative Dq. 

The terms of first order in are introduced exactly as in [1], but we use a 

slightly different notation which is convenient when one does not go beyond that order. 
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2.1 Formulation 



The lowest order (static) Lagrangian, 

£ s tatO) = iphix) A) ih(x) , (2.1) 

is written in terms of the backward covariant derivative Dq as in [27] and the 4- 
component heavy quark field subject to the constraints P+iph = V'h > V'h-ff = V'h W1 ^h 
P + = (1 + 7o)/2. At the first order we write the HQET Lagrangian 

£hqet(x) = £ s tat(x) + £ (1) (x), (2.2) 
C^\x) = -uj spin O spin (x) - uj kiQ O kiQ (x) , (2.3) 
Cspin(x) = ^ h (x)vB^ h (x) O k in = 4WDVhW, (2.4) 

such that the classical values for the coefficients are Wkm = w sp m = l/(2m\ ) ). We use 
the discretized version er-B = Ylk j a kjFkj/(2i) , with cifcj and the lattice field tensor F 
defined in [28]. The kinetic term D 2 is represented by the nearest neighbor covariant 
3-d Laplacian. The effective theory is expected to be renormalizable at each (fixed) 
order in l/my, if (and only if) path integral expectation values are defined by expanding 
the path integral weight as [1] 

exp(-q 4 y^[£ H QET(^) + Aight(ff)]) = exp(-a 4 ^[£ stat (x)+£i ight (x)]) (2.5) 

X 

1 + a 4 ^2[u sp mO spin (x) + WkinOkin(^)]) • 



x 



For correlation functions of some multilocal fields O this means 

(O) = (0) s tat + Wkm« ^(OOspin(x))stat (2.6) 

x x 

= (C)stat + Wkin(O) kin +^spin(C ) spin , (2.7) 

where (0) stat denotes the static expectation value with Lagrangian Atat(aO + Aight 
All terms composed of just the relativistic quarks and the gauge fields are summarized 
in Aight (a:)- Note that as one performs the Wick contractions of the heavy quark field, 
the 1/mb terms O k \ n (x) } O sp i n (x) leave behind insertions in the static heavy quark prop- 
agators. From the point of view of renormalization all terms in eq. (|2.6p are correlation 
functions in the static effective theory, which is power counting renormalizable. 

The above form assumes that O contains all 1/rrtb terms needed to represent the 
local fields in the effective theory. A relevant example is the time component of the 
heavy light axial current. In the effective theory it is represented as 

Aq QET (x) = ZA QET [^ s o tat (x) + CA QET M s o tat (x)], (2.8) 
A^\x) = ^ x {x) l0l ^ h {x) , (2.9) 

5A^\x) = ^(^(Vi+Vfh^h^)- (2.10) 
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Later we will also use the space components of the vector current represented by 

Vf QET (x) = Z^ QET [F fc stat (x) + c^ QET ^ stat (x)], (2.11) 
Vt\x) = ^(x) 7 ^h(x), (2.12) 

5V^\x) = -^(x)i(V i+ V*) 7j7 ^h(x). (2.13) 

We have chosen V fc stat , 5V fc stat such that they are exactly related to Aq^, by a spin 

rotation. 

The coefficients u>ki n , w sp i n , Z^ ET , c^ ET , Zy I< ^ ET , Cy^ ET are functions of the bare 
coupling g$ and of the heavy quark mass in lattice units. They represent bare parameters 
of the effective theory, which are to be fixed by matching to QCD. Just like u;ki n! w sp i n , 
the coefficients c A , c v are of order l/m^,, while we may write 

Z HQET = ^stat + ^(.1) ^ tfl) = (l/ mb ) , ( 2 .14) 

and similarly for Zy^ ET Q. Note that in the expansion to first order, terms such as 

^kin c A ^ V m b are *° ^ e dropped. 

Below we will consider an example and discuss that indeed the bare parameters 
Wkm ; ^spin 5 ^A^ ET ' c A T '^v T ' C V^ ET ana - m bare are sufficient to absorb all diver- 
gences in the effective theory at this order in 

2.2 1/wib expansion in a geometry without boundaries 

In order to illustrate further how the expansion works, we consider a two-point function 
of a composite field in a space-time without boundaries, i.e. with periodic boundary 
conditions or in infinite volume. We choose the example 

C AA (xo) = Z 2 A a 3 £ (a (x)(A )\0)} (2.15) 

X 

with the heavy-light axial current in QCD, A^ = V , i7/^75'0b) an d Z A ensuring the nat- 
ural normalization of the current consistent with current algebra [29,30]. The l/m^ 
expansion reads 

Caa(^o) = e- mb -^(Z™) 2 a 3 ^[(^ tat (x)(^ tat (0)) t )stat (2.16) 

X 

+ u; kin (^ tat (z)(^ tat (0))t) kin + a; spin (^ s tat (x)(^ tat (0)) t )s P in 

+ C^ QET (^ S tat (x)(Mr t (0)) t ) s tat + C? QET (5< at (x)(^ tat (0)) t )stat 

= e- m "(^A QET ) 2 [cfi*(xo) + u kin C^(x ) + u^Clfixo) 

+c» QET [Clil(x ) + CtH(xo)]\ (2.17) 



1 If O(a) improvement is desired in the static approximation, there are also aSA^^, a<5V fc stat cor- 
rections to the currents. They are not relevant in the present discussion but will be taken into account 
when necessary. 
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up to terms of order As mentioned in the introduction, the mass shift mbare = 

O(m-b) includes an additive mass renormalization. It is also split up as 

mbare = Kare + ™bare > with ™bare = 0(l/m b ) , (2.18) 

and the expansion e -mbaro:ro = e~ m ^^ x °(l — ^o m bare) * s understood. 

For illustration we check the self consistency of eq. (|2.17p . The relevant question 
concerns renormalization, namely: are the "free" parameters m-bare • • • c^ ET sufficient 
to absorb all divergences on the r.h.s.? We consider the most difficult term, C^(xq). 
According to the standard rules of renormalization of composite operators, it is renor- 
malized as 

fe) (so) = e-^^^fa^^W (^ tat (0)) f (O kiR ) R (z)) + C.T. , 

\ f JrL » ' Stilt 

X, z 

(2.19) 

where C.T. denotes contact terms to be discussed shortly. The renormalized operator 
(Okin) R (-z) involves a subtraction of lower dimensional ones, 

{O kin ) R (z) = Z 0kin (O kin (z) + ^ h (z)D ?J; h (z) + ^ h {z)^{z)) , (2.20) 

written here in terms of dimensionless c%. Since we are interested in on-shell observables 
(xq > in eq.(2.19)), we may use the equation of motion DqiP^(z) = to eliminate the 
second term. The third one, ^^^(z)^h(z), is equivalent to a mass shift and only changes 
m bare> w hich is hence quadratically divergent El Thus all terms which are needed for the 
renormalization of Okin are present in eq. (|2.17p . It remains to consider the contact terms 
in eq. (|2.19p . They originate from singularities in the operator products O k i Q (z)A ! Q at (x) 
as z — > x (and O kiri (z) (^4Q tat )^(0) as z — > 0). Using the operator product expansion 
they can be represented as linear combinations of A^ a,t (x) and <5AQ tat (x). Such terms 
are contained in eq. (|2.17|) in the form of C^jf , Cf^ and C^f^ H 

We conclude that all terms which are needed for the renormalization of C^™(xo) 
are present in eq. (|2.17p ; the parameters may thus be adjusted to absorb all infinities 
and with properly chosen coefficients the continuum limit of the r.h.s. is expected to 
exist. The basic assumption of the effective field theory is that once the finite parts of 
the coefficients have been determined by matching a set of observables to QCD, these 
coefficients are applicable to any other observables. 

2 Using the explicit form of the static propagator, eq. (2.4) of reference [27], one can check that indeed 
a 3 E x (^o tat (z) (^o tat (0)) t a 4 E z ?h(^)V'h( 2 )) = XoCff(x ). 

3 t t stat 

Aq a (a;) and SAq a (x) are the only operators of dimension 3 and 4 with the correct quantum 
numbers. Higher dimensional operators contribute only terms of order a. Note that the Aq (x) term is 
power divergent ~ l/{am- b ). This divergence is absorbed by a power divergent contribution to Z^ ET 
(at order 1/w.b). 
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The B-meson mass is given by Caa(^o) i n large volume via 

m B = - lim ^±^0 logCAA(xo), (2.21) 

XQ-^fOO Z 



with 



dof(x ) = -[f(x + a) - f(x )] , d* f(x ) = -\f(x ) - f(x - a)} . (2.22) 
a a 



Inserting the HQET expansion we derive 

m B = m g at + mjj' , (2.23) 

with 

< at = m^ e + E Btat , E stat = - lim log (*„) , (2.24) 

xo— >oo Z 

m$ = + ^kin-Erin + WspinE'spin , (2.25) 

E kin = - lim ^±^- [C%Z(x )/Cft(xo)] , (2.26) 

E spin = - x ^^^[C^(x )/CtX\x )]. (2.27) 

Here the terms oc c^ ET of eq. (|2.17|) do not contribute. They are proportional to 
the derivative of ratios C| A a A (xo)/C A t A *(xo). At large xq these ratios approach a con- 
stant since SA^^ has the same quantum numbers as j4Q tat . Using the transfer matrix 
formalism (with normalization (B\B) = 1), one further observes that 

Ekin = -(B\a? Ckm(0, z)\B) stat , E spin = -(£|a 3 ^ O spin (0, z)^}^ . (2.28) 

z z 

As expected, only the parameters of the action are relevant in the expansion of a hadron 



mass. In the above relations absorbs a linear divergence of E stat and rn^ re 



quadratic divergence of E kiD . 

Going through the same steps in the vector channel and using the spin symmetry 
of the static action is one way to see that the combination 

m B = \\. m ^ + 3m B*] = Wbare + E sta t + W kin E kin (2.29) 

is independent of uJ S pin- It is instructive to represent this equation in a different way, 
subtracting the 1/a (and I /a 2 ) divergences of E stat (and -Ekin)- In this way we have 

mf = mr'+mP+mP + mf , (2.30) 



fOo) 



+ , (2.31) 



= E stat -E b s t u £, (2.32) 
< a) = m^+u^E^, (2.33) 
4" = ^kin[E kin -E^ n b ], (2.34) 
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with finite terms rrtg ,mn ,mj ,mg . Our strategy, described in the introduction 
can be seen as a way of determining the coefficient as well as the subtractions 
i^ltat, -^Sn fr° m hriite volume computations in QCD and HQET. Finite parts in the 
subtraction terms do of course depend on the detailed choice of kinematical parameters 
such as the matching volume, but the end result is unique up to terms of order 1/m?. 
Note that by the same logics, the order 1/mb term, mg 1 "' + m^ b \ is not unique but 
depends on the details of the strategy. 

Since the prediction eq. (12,29h requires only the knowledge of two parameters, we 
also need only two finite volume observables to perform the matching with QCD. The 
Schrodinger functional is particularly useful to find suitable observables [1,22,31]. We 
proceed to discuss the 1/mb expansion in this situation. 



2.3 Schrodinger functional correlation functions 

The pure gauge Schrodinger functional has thoroughly been discussed in [24], relativistic 
and static quarks were introduced in [25] and [28]. In particular in the last reference 
also Symanzik 0(a)-improvement was discussed. The improvement of the Schrodinger 
functional requires the addition of dimension four local composite fields localized at 
the boundaries [32]. However, it turns out that there are no dimension four composite 
fields which involve static quarks fields and which are compatible with the symmetries 
of the static action and the Schrodinger functional boundary conditions and which do 
not vanish by the equations of motion. Thus there are no 0(a) boundary counter terms 
with static quark fields. For the same reason there are also no 0(l/mt,) boundary terms 
in HQET. This then means the HQET expansion of the boundary quark fields £, C is 
trivial up to and including terms. 

For details of the boundary conditions as well as the definition of the fields £ 
we refer to [28], where also our notation is explained. For a general understanding it 
is, however, sufficient to note a few facts. In space the fermion fields are taken to be 
periodic up to a phase, 

ip(x + kL) = e i6 ip(x) , ^(x + kL) = e^V^) , (2.35) 

with the same phase 9 for all quark fields, whether relativistic or described by HQET. 
In time we take homogeneous Dirichlet boundary conditions at xq = and xq = T [28]. 
Correlation functions can be formed from composite fields in the bulk, < xq < T, 
and boundary quark fields C > C- I n QCD, correlation functions in the pseudoscalar and 
vector channel are 

fA(xo,0) = -y^((Ai) ( ; r)Cb(y)75Ci(z)) J (2.36) 
k v (x 0: 9) = ^ <(y I ) fc (x)C b (y) 7fc Ci(z)> • (2.37) 
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The 0(a) improved currents Ax, V\ can be found in [1]. Furthermore we consider bound- 
ary to boundary correlation functions 

tm = E <Ci / (u) 7 5C'b(v)C b (y)75Ci( Z )> , (2.38) 

u,v,y,z 

ki(9) = ~ Y, (Cl'CuhfcCUvKbtyb^)) - (2.39) 

u,v,y,z,fc 

Their renormalization is standard [33], for example 

[/a] r (x , 9) = Z A Z 2 c f A (x , 6) , [fi] R (9) = Z\h{9) , (2.40) 

with Z^ a renormalization factor of the relativistic boundary quark fields. 

In complete analogy to the case of a manifold without boundary we can write down 
the expansions to first order in 1/mt,. They read 

[/a] R = ^ QET ^Ce~ m ^ ,(2-41) 
[ky] R = Z^ ET Z Ch Z c e- m ^° )(2 . 42) 

-^ QET Z Ch Z C e- m <— o jjstat + c HQET jstat + ^jUn _ 1^/^*} , 

[/i]r = ZlZle- m ^ T {f^ t + u^f^ + u spin ff n } , (2.43) 
[h] R = ^Zle-^j/r+^/^-i^/r 11 } • (2.44) 
Apart from 



/fi at (x o ,0) = -f E <^o tat (^)C h (y)75Ci(z)> (2.45) 



2 



y 5 z 



the labeling of the different terms follows directly the one introduced in eq. (|2.7p . We 
have used identities such as /^ in = — ky m , /^ pm = 3ky in . As a simple consequence of 
the spin symmetry of the static action, these are valid at any lattice spacing. 

3 Finite volume observables and step scaling functions 
3.1 Observables 

We concentrate on a strategy based on the correlation functions fx,k\ alone. This is 
advantageous, since the additional coefficients c E( ^ ET , Cy^ ET in eq. (|2.41j) . eq. (|2.42p are 
avoided. Apart from the b-quark, we set the masses of all quarks to zero. 
In terms of the spin-averaged combination, 

F 1 (L,9) = i[log/i(0) +31ogfci(0)] , (3-1) 
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we form 



R 1 {L,e 1 ,9 2 ) = F 1 (L,8 1 )-F 1 (L,8 2 ) at T = L/2 (3.2) 
r^L,^) = - dT ^ d * T F 1 (L,9 )atT = L/2. (3.3) 

Note that the boundary quark wave function renormalization cancels in Ri and in T\. 
They are thus finite after renormalization of the parameters of the Lagrangian. 
The dimensionless observables, 

*i(L,M b ) = R 1 (L,e 1 ,e 2 )-Rf^(L,e 1 ,e 2 ), (3.4) 

$ 2 (L,M b ) = LriGMo), (3-5) 
Rf*\L,e h e 2 ) = log [fl ta,t (L,9i)/ fl ta,t (L,8 2 )] atT = L/2 (3.6) 

are parametrized in terms of the RGI mass of the b-quark, Mb. They have a particularly 
simple 1/fTib expansion 

*i(L,M b ) 
c& 2 (L,M b ) 

which involves 

R\ m (L, 9i,9 2 ) 

rf at (L,9 ) 
rf n (L,9 ) 

The 9q,9\,9 2 dependence of 3>i is not explicitly written, but will of course be relevant 
in the numerical results. For the reader familiar with [1,34], we point out that I\ 
differs from T which was used in those references. Note that in eq. (|3.4p we subtract 
the static term. This simplifies subsequent formulae. In fact, whenever such a lowest 
order contribution is universal (in the sense of having a universal continuum limit) and 
independent of an HQET parameter, it will be convenient to subtract it. Despite this 
subtraction, we refer to $1 as an observable in QCD. 

The reader may be surprised that we introduce the quantity T± which contains a 
(discretized) derivative with respect to the time extent, T. Its MC evaluation requires 
two separate simulations @. However, obviously a quantity of order is needed and 
this is obtained from some logarithmic derivative of a correlation function. Boundary- 
to-boundary correlation functions are then very convenient since one does not have to 
deal with the l/m^ corrections to the currents. It is a useful feature of the Schrodinger 
functional that such gauge invariant correlation functions are available. 

4 In App.[C]we discuss a different strategy, which is based on the xo-derivative of /a and thus requires 
less simulations. Note, however, that these additional simulations do not represent a significant effort. 



UkmRi m (L,9i,9 2 ) , (3.7) 
L [m bare + Tf at (L, 8 ) + Wkinlf n (L, 8 )] , (3.8) 



fi in (L, 9\) ff D (L,9 2 ) 



atT = L/2, (3.9) 



fl tat (L,9t) ft t&t {L,8 2 ) 

logjf^o) atT = L/2, (3.10) 

■ltf m (9 )/f! tat (8 )} atT = L/2. (3.11) 



8 T + <)j r: M 



9 



3.2 Step scaling functions 



We turn to the relations between $j(L,Mb) and <5j(2L,Mb) in the effective theory. 
The dimensionful variable L is replaced by the Schrodinger functional renormalized 
coupling g 2 (L) [35] over which we have good control in numerical computations [20]. 
Straightforward substitution yields 



$!(2L,M b ) 
$ 2 (2L,M b ) 



af n (u)^(L,M h ), 

2$ 2 (L, M b ) + a m (u) + af a {u) $ X (L, M b ) 



(3.12) 
(3.13) 



where always u = g 2 {L). Our continuum step scaling functions a (with any subscripts 
or superscripts) are defined in terms of those at finite lattice spacing as 



a(u) 



lim S(u, a/L) 

a/L^O 



At finite lattice spacing we have 



u, a/L) 



R\ m (2L, 61,62) 



R 1 n (L, 61,62) 



u=gl(L) 



^T(u,a/L) 
H m (u,a/L) 



2L [T^(2L,6o)-T^(L,6 )} 



Rf»(L,9 1 ,e 2 ) 

l *(2L,0 o )-rf 



2L [Tf- t (2L,6 )-Tf-\L,6 )] u 



u=g 2 {L) 
9 2 (L) 



(3.14) 

(3.15) 

(3.16) 
(3.17) 



The above equations are easily derived. In a first step, just from the 1/my, expansions of 
one obtains them at a given resolution a/L or equivalently at fixed bare coupling, g$ . 
One then uses that &i(L, Mb) are dimensionless physical observables with a continuum 
limit. Since the self energy of a static quark cancels in a m , also that quantity has a 
finite continuum limit. Thus the continuum limit of the step scaling functions S m , 
exists and eqs. (|3.12l3.13p can be written in terms of continuum quantities, as we have 
done. 



4 M b including 1/mb corrections 

Before giving the equation for Mb, we recall the overall strategy. For L\ ~ 0.4 fm 
we compute <3?i(Li, Mb), <J>2(-^i, Mb) for a few quark masses around the physical one in 
quenched QCD. It is understood that the continuum limit is reached by an extrapolation 
and with a suitable interpolation of <!>j in Mb, these quantities can be considered to be 
known as a function of Mb. With the step scaling functions described in the previous 
section and computed in the effective theory, we then arrive at $1(^2, Mb), $2(^21 Mb), 
where L2 = 2L\. It remains to express the spin averaged B-meson mass m| v in terms 
of $ 1 (L 2 ,M b ) J $ 2 (L 2 ,M b ). 
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To this end, we straightforwardly combine eqs. (|3.7|3.8p with eq. (|2.29j) and obtain 

Lmf = $ 2 (L, M b ) + L[£ stat - rf '(i, 9 )] + ^—1^$^, M b ) , 

(4.1) 

We now set L = L 2 in this equation and insert eq. (|3. 13j) . In the form of eq. (|2.30p we 
then have 

L 2 m ( ° a) (M h ) = a m (n 1 ) + 2$ 2 (L 1 ,M b ) (4.2) 
L 2 mf ] = L 2 [E st ^-Tf^(L 2 ,0 o )], (4.3) 
L 2 m^ a) {M h ) = Sx^Mb), (4.4) 

pkin "pkin It a \ 

L 2 m^ b \M h ) = L 2 L l ™ \ 0) of (m) M h ) , (4.5) 

-Kl (^2 > "l j "2 J 

where 

ui=5 2 (^i), L 2 = 2Ll (4.6) 



The subtraction of power divergences in eq. (|2.32j) . eq. (|2.34j) are £^ = rf at (L 2 ,0o), 
= Tf n (L 2 ,e ) and of n («i) M h )/R\ m {L 2 , 9 U 2 ) is a representation of the 

bare parameter u;ki n in eq. (|2.34[) . The other parts, mg"',™^, are computable entirely 
in finite volume. 

The step scaling functions a have been discussed before. They can be computed 
with lattice spacings such that a/L\ is reasonably small, say below 1/6. Of course they 
should be extrapolated to the continuum. We work with lattice spacings a < 0.07 fm 
in this step. The relativistic observables $j(Li,Mb), i = 1, 2 are computed for a < 
0.02 fm, where a relativistic b-quark can be described by the 0(a)-improved Wilson 
action with controlled a 2 -effects. Finally, the combinations L 2 [E 8ta,t — Tf at (L 2 , 6q)] and 

L 2 nkin ,r „ ' , are computed for lattice spacings of a < 0.1 fm such that finite size 
effects in E stat and E^ n are negligible on lattices with an affordable number of points. 
The mass of the b-quark is obtained from eq. (|2.30[) by expanding 

M b = M^ 0) + , (4.7) 
where Ml is the solution of the static equation 

m - = r4° a) {M® ) + r4 0b) {M® ) (4.8) 
and the 1/m-b correction is 

M« = -^4 la) K 0) )+4 lb) K 0) )] (4-9) 
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Figure 1: Continuum extrapolation of $2(^1, Mh), for z = 10.4 , 12.1 , 13.3 from bottom to top. The 
errors in the relation between bare quark mass m q and the RGI mass M are translated into errors in 
$2- The go-independent part of that error is included after [31] the continuum extrapolation (left side 
error bar). On the right, the equivalent in the alternative strategy is shown for 0q = 1/2 (see App. [Cj. 



with 



S = ^[mf)(M b )+mf)(M b )]=^-[4 0a )(M b )]. (4.10) 

We finish the discussion of the strategy with a remark on the dependence on the 
mass of the light quarks. This is relevant because it is of course better to consider the 
spin-averaged B s quark mass in eq. (|2,29p ; the necessary large volume computations are 
easier than for the meson. In the quenched approximation the parameters in the ac- 
tion rribare, <^kin are independent of the light quark massH Since our strategy determines 
them through finite volume computations, it follows that in all these computations the 
light quark mass may be set to zero, a convenient choice. Only and -E s tat are then 
to be computed at the mass of the light quark of the meson who's (spin averaged) mass 
is considered. 

5 Results 

We have performed a numerical computation in the quenched approximation, using the 
O(a) improved Wilson action [32,41,42]. The box size Li is chosen as L2 = 1.436ro, 



5 In general, 5m (and hence also mbarc) will contain a term like b(go)m\, where for simplicity the light 
quarks are assumed to be degenerate with mass m\. Obviously, b(go) — O(po) does, however, vanish for 
N{ = 0. As a renormalization term odd in mi, it is also absent for twisted mass lattice QCD [36] and 
QCD with exact chiral symmetry [37-40]. 
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(a/L,) 2 (a/Lj) 2 

Figure 2: Continuum extrapolation of $i(Li,Mb), separately for R\(L%, 1/2, 1) in QCD (left) and 
for J?f at (Li, 1/2, 1) in the static approximation (right). Circles denote results with action HYP1 and 
squares, displaced slightly for visibility, are from action HYP2. The corresponding continuum extrapo- 
lation lines are slightly displaced as well. 

where ro, defined in terms of the static quark potential [43], has a phenomenological 
value of ro ~ 0.5 fm. From [20] we know the Schrodinger functional coupling g 2 (L\) = 
g 2 (L2/2) ~ 3.48. Given the knowledge of r^/a as a function of go of Ref. [44] and that 
of the renormalized coupling [20], it is then convenient to fix go in different ways for 
the different steps of the calculation. The differences are of course only a-effects which 
disappear in the continuum extrapolations. We give more details below. We will take 
the uncertainties in the relations g 2 (Li) ~ 3.48 and g 2 (Li/4) ~ 1.8811 (which we need 
later) into account in the very end. 

In order to complete our definitions, we further choose 8q = and 9±, 62 £ {0, 1/2, 1}. 
The different values of 9%, 82 offer the possibility to check whether our final results are 
independent of these arbitrary parameters as they should be up to small l/m 2 , terms. 
Simulation parameters as well as raw results are listed in tables in App. |A] and [Bj 

5.1 QCD observables 

For this part of the computation, we determined the bare parameters as in [31]: go is 
fixed by requiring g 2 (L\/4) = 1.8811 for given resolutions a/L. The PCAC mass of the 
light quark, defined exactly as in that reference, is set to zero. Our heavy quark masses 
are chosen such that z = L\ ~ 10 — 13 . The bare parameters are listed in Table [3l 
We focus our attention directly on the continuum extrapolations. As an example 
we show $i(M b ,Li) and ® 2 {M h ,L 1 ) in Fig. [2] and Fig. CD Note that for the static 
subtraction i?^ tat (Li, 1/2, 1), displayed on the right of Fig. [21 our lattice spacings are 
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roughly a factor three larger, since in the effective theory we only have to respect 
a/L\ <C 1, not aMb *C 1 (for details see App. [A]). Data have been obtained for two 
static actions, HYP1 and HYP2 [27]. In fitting them to the expected a-dependence, 
their continuum limit value is constrained to be independent of the action, but the a 2 
slopes are of course different. The data for the different actions are highly correlated. As 
in all such cases, the errors of the continuum limit are computed from jacknife samples. 

For values of 6i which differ from the choice made in the figures, the a-dependence 
is very similar. In all these cases we find that extrapolations linear in a 2 using all four 
available lattice spacings are compatible with the ones where the data point at largest 
lattice spacing is ignored. We take the extrapolations with three points as our results 
for further analysis, since they have the more conservative error bars. The continuum 
limits are listed together with the raw numbers in Tables [6] and [H From a fit of the 
continuum ^(-z) to a linear function, we then extract the slope 

S = = o.61(5) (5.1) 

and we are done with the matching. The rest of the numerical computations is carried 
out in the effective theory. 

5.2 HQET step scaling functions 

Next we discuss the connection of §>i(Li, M^) to $>i(L 2 , M^), L 2 = 2L\. It is given by 
the step scaling functions of Sect. I3.2L The bare parameters used in their computation 
are described in App. [Aj, and the values at finite resolution a/L\ are given in Tables 

[Hoi. 

At lowest order in l/mt>, only cr m contributes. In its continuum extrapolation 
(Fig. El Table [8]) we allow for a slope in a 2 , although the data are compatible with a 
vanishing slope. Note that the absolute error of <r m is negligible in comparison to twice 
the one of $2 (see Fig. [I]) to which it is added in eq. (|4.2p . In fact the uncertainty in 
cr m corresponds to an error of only 5 MeV in the b-quark mass, illustrating the possible 
precision in the static effective theory with these actions [27, 45] . 

A relevant question is how the precision deteriorates when one includes the first 
order corrections in l/mt>. Then two more step scaling functions contribute. In Fig. HI 
we illustrate how the continuum limit of a\ m is obtained. Here we have to allow for a 
linear dependence on the lattice spacing, since the theory is not 0(a) improved at the 
level of the contributions [1]. Taking the more conservative fit with only three 

points, we arrive at the continuum limit listed in Table [9] for all combinations #i,#2- 
In eq. (|4.5I) . a\ m is multiplied by small numbers (of order 1/mb). This means that its 
error will be negligible in the overall error budget. 

Instead of a^™ we show directly the continuum extrapolation of m^ a \ eq. (14, 4h . 
As for E m , the data shows no significant a-dependence. Nevertheless, in order to have 

6 For Em, S2" 1 the coarsest resolution considered is a/Li = 8. Due to the derivative dr at T = L/2, 
smaller values of Li/a would involve a very short time separation. 
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(a/Lj) 2 

Figure 3: Continuum extrapolation of E m and E : 



a realistic error estimate, we allow for a linear slope in a (Fig. [5]). In Table [TOl we list 
the raw numbers for S| m as well as the extracted continuum limit for further analysis. 

5.3 Large volume matrix elements and Mb 

The last missing pieces in eq. (I2.30|) are the large volume static energy E s tnt , eq. (I2.24J) , 
and the matrix element of the kinetic operator E^[ n , eq. (|2.28p . Here, in contrast to the 
rest of our numerical evaluations, the light quark mass is set to the mass of the strange 
quark in order to avoid a chiral extrapolation. The spin averaged mass of the B s system 
is then to be inserted into eq. (|2.29p . 

Although -E s tat and -Ekin can be computed with periodic boundary conditions we 
here follow [46] and evaluate also these quantities with Schrodinger functional boundary 
conditions in a large volume of about T x (1.5fm) 3 , with 1.5fm < T < 3fm (also a check 
for finite size effects is carried out). The extraction is fairly standard, but still care 
has to be taken to make sure that the ground state contribution is obtained. This is 
a particularly relevant issue for B-physics, because the gap to the first excited state 
is rather small. We relegate details to App. [B] and discuss immediately the universal 
combinations L2m^ b \ L2rn^ which enter in eq. (|2.30p . The static contribution, shown 
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a/Lj a/Lj 

Figure 5: Lattice spacing dependence of L^rn^ 1 for Mb = M^°\ On the left we show mg°' as 
introduced in Sect. [4] with 6\ = 1/2, 6*2 = 1. We insert $i in the continuum limit, such that the 
lattice spacing dependence is just due to E3 • O n the right the corresponding quantity is shown for 
the alternative strategy of App. [C] again with continuum values for $i. There we set 6a = 1/2 , 9\ — 
1/2 ,0 2 = 1. 
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in Fig. [6l is known with very good precision |j. In contrast, the l/rrtb correction L2 m B 
does have a noticeable total uncertainty (Fig. El Table [Q). Still, this error is only about 
50% of the one on 2<&2- Note also that this error is almost entirely due to E^[ n which 
may possibly be computed more precisely by other techniques [47]. 

We now have all pieces necessary to solve the equations for Mb. The static one, 
eq. (|4.8f) , is illustrated in Fig. [HJ The horizontal error band is given by subtracting the 
static pieces a m + i^^R from the experimental number 

m | v = 5.405 GeV . (5.2) 

The figure demonstrates again that the main source of error is contained in the QCD 
computation of <3?2- Finally, by interpolating $j(Li,Mb) to Mb = M^ we obtain 
(#i = l/2,0 2 = l) 

17.25(20) , M^ 0) = 6.806(79) GeV for r = 0.5 fm (5.3) 
-0.12(9) , M^> = -0.049(39) GeV for r = 0.5 fm (5.4) 
17.12(22) , M b = 6.758(86) GeV for r = 0.5 fm . (5.5) 

Here the small difference <? 2 (Li/4) — 1.8811 as well as the statistical uncertainties in 
g 2 (L\) and L\/tq have been taken into account, as explained in App. [Dj Moreover, 
one can see in Table [1] that the 9i dependence of the l/m-b contribution is absorbed. 
With Ajyjgro = 0.602(48) [35,48], the 4-loop (3 function and the mass anomalous dimen- 
sion [49-52], we translate Mb = + to the mass in the MS scheme at its own 

7 We show the results given for the static action HYP2. The continuum extrapolation with action 
HYP1 looks very similar, but the fit has a smaller \ 2 - 
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a/L 2 a/L 2 
Figure 7: Lattice spacing dependence of L2TOg , details as in Fig. [5] 
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1/2 


-0.06(3) 


-0.06(8) 


1/2 


1 


-0.06(3) 


-0.06(8) 


1 





-0.06(3) 


-0.06(8) 



Table 1: RGI results of 1/mb correction of the b-quark mass, in units of ro. 

scale, 

fn h {fn h ) = 4.347(48)GeV ; (5.6) 

the associated perturbative uncertainty can safely be neglected. For completeness we 
note that in the MS scheme the 1/mb term amounts to — 27(22)MeV. 

5.4 Comparison to results from an alternative strategy 

As mentioned earlier, at first sight it appears more natural to base the computation 
of Mb on the logarithmic derivative of the spin average of /a and ky as the prime 
finite volume quantity. We have not chosen this option as our standard strategy since 
then three observables are needed for matching. However, it is useful to consider also 
that alternative strategy in order to perform an explicit check that l/m-j^ terms are as 
small as expected. The results can be appreciated without detailed definitions of the 
observables and step scaling functions, the interested reader can find them in App. [Cj 
Here we note that within this alternative strategy we actually have nine different sets 
of {$i, <3?2; ^3}- Only one observable, <J>i = <E>i, is in common to the two strategies. For 
our graphs we have selected (arbitrarily) one choice of parameters. 
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0o = O = 1/2 O = 1 





1/2 


17.12(22) 


17.25(28) 17.23(28) 17.17(32) 


1/2 


1 


17.12(22) 


17.23(27) 17.21(27) 17.14(30) 


1 





17.12(22) 


17.24(27) 17.22(28) 17.15(30) 



Table 2: RGI results of Mb inlcuding the \/mh correction, and comparison of the two strategies. 



First, let us summarize what kind of differences one expects in such a comparison 
apart from a-effects. In the order of magnitude counting, we take L^ 1 ~ A ~ 0.5 GeV 
and of course = 2L\. The matching observables Ti, T av are constructed to be equal to 
the quark mass at the leading order in the HQET expansion. They start to differ at the 
next to leading order, which means by terms of order A. Also their dependence on 8i is 
of that magnitude. Since ^(-^li Ma) and $3(^1, Mb) have been made dimensionless by 
multiplication with L\ and L\ happens to be around A -1 , the differences of $2(^1, Mb) 
and ^(Li, M3) are order one. The step scaling functions cr m , <r m as well as L 2 mg*> arc 
added to $2 (or $3) to obtain I^me i n static approximation. Thus they depend on 
the details at the same level, apart from a trivial LijL\ = 2 factor. Of course, in the 
total static estimate tqM^ this dependence is reduced to r^M^ x (A/mb) 2 ~ 1/5. In 
the same way, the l/rn^ corrections L2m!^ a \ I^rn^ themselves have a dependence on 
the matching conditions which is L2 x A 2 /m^ ~ 1/5 but the final result r^M^ including 
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these terms is accurate and unique up to ro x (A 3 /mt,) 2 ~ 1/50 corrections. 

We leave it to the reader to check in Fig. [1] to [7] that these expectations are fully 
satisfied by our results El- In fact it appears that our estimate for the expansion pa- 
rameter, A/mb ~ 1/10 is quite realistic. Of course, to find this out requires an explicit 
computation of the correction terms as presented here. In some cases, such as rag 1 ''', 
our precision is not good enough to resolve a dependence on the matching conditions. 

In the b-quark mass in the static approximation, roM^ (eq. (|5.3f) and Table fT2|) . 
the maximum difference is 0.5(2), which is of the predicted order of magnitude. Finally, 
when we add all contributions together, the results from the alternative strategy, Tabled 
are fully in agreement with eq. (15.5f) . As expected 1/m 2 terms are not visible with our 
precision. They can safely be neglected. 

6 Conclusions 

The main conclusion of this work is that fully non-perturbative computations in lattice 
HQET, as they have been suggested in [1], are possible in practice. In particular, the 
uncertainties in the l/m,b corrections are smaller than those in the static approximation, 
despite the fact that we numerically cancel large a~ 2 divergences in the 1/mb terms. 
The final error in the mass of the b-quark is dominated by the uncertainty in the 
renormalization in QCD. Errors due to simulations in the effective theory can almost 
be neglected in comparison. 

A very nice result is the independence of the final numbers for Mb of the matching 
condition: Table [2] shows that within our reasonably small uncertainties, we get the 
same results for the quark mass for altogether twelve different matching conditions. 
This is expected up to very small terms of order tqM^ x (A/mb) 3 ~ 0.02, which should 
be compared to our result tqM^ = 17.12(22) — 17.25(28). Here the quoted range is due 
to the different matching conditions. In the order of magnitude estimates we have made 
a guess for the typical scale of A ~ 0.5 GeV. In the static approximation, some of the 
matching conditions yield slightly differing results for the quark mass in agreement with 
the expectation for such variations of tqM\> x (A/mb) 2 ~ 0.2. 

Both this explicit test of the magnitude of the different orders in the expansion 
and the naive order of magnitude estimate say that 1/m^ corrections are completely 
negligible. 

Still, in aspects of the computation, considerable improvement can be envisaged. 
For example, return to the 1/mb contribution to the B-meson mass Fig. The statis- 
tical errors grow rapidly as one decreases the lattice spacing. The by far dominating 
uncertainty in the shown combination is the one of the large volume matrix E^in- It ap- 
pears worth while to look for improvements, maybe along the line of [53]. Due to these 
errors, and of course the missing 0(a)-improvement of the theory at order 1/mb [1]> 

8 We note in passing that E m , in contrast to E m , does in principle require an improvement coefficient, 
c\ at [28], for 0(a)-improvement. It has been set to the 1-loop values from [27], but the results are rather 
insensitive to cj at , so its uncertainty can be neglected. 
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the continuum extrapolation is not easy. Fortunately it is still precise enough for the 
present case. It will be very interesting to see cases where the l/m^ corrections are 
larger, as it is expected, for example, for 

Let us now turn to the computed value of Mb, eq. (15. 6p . Starting from a precisely 
specified input, namely ro, tok and (mB s + 3mB*)/4, the value of M b is unambiguous 
in the quenched approximation, because these inputs fix the bare coupling, strange 
and beauty quark masses. We have used the experimental meson masses and ro = 
0.5 fm. Our numbers for M b or fn^ may then be used as a benchmark result for other 
methods. Indeed, a comparison shows agreement with [54] and the recent extension of 
that work [55] m b = 4.42(7) GeV. 

Earlier, the review [56] quoted fn h = 4.30(5)(5) GeV and m b = 4.34(3)(6) GeV, 
based on static computations [57] and an extrapolation of NRQCD results to the static 
limit [58] respectively. A perturbative subtraction [18, 59, 60] of the linear divergence 
5m was carried out in these static estimates and, of course, a continuum extrapolation 
could not be done. 

However, if other inputs are used, the result may change because ro is only ap- 
proximately known and because the quenched approximation is not real QCD. A rough 
idea on the possible changes can be obtained by varying ro by ±0.05 fm. This changes 
^b(^b) by roughly ±80MeV. 

These remarks just serve to stress the obvious necessity of performing computations 
with Nf > 0. The ALPHA-collaboration is presently starting with Nf = 2, where 
the renormalization of the quark mass in QCD is known [21]. The necessary HQET 
computations are not expected to be a big numerical challenge, apart from the large 
volume B- meson matrix elements: simulations of the Schrodinger functional for L < 1 fm 
are not very demanding with nowadays computing capabilities [61]. Altogether the 
extension of the present work to full QCD is feasible and should be carried out, since 
presently no better method is known to compute the b-quark mass from lattice QCD. 

Acknowledgements. We thank Stephan Diirr for collaboration in the early stages 
of this work [62]. We thank NIC for allocating computer time on the APEmille com- 
puters at DESY Zeuthen to this project and the APE group for its help. This work is 
supported by the Deutsche Forschungsgemeinschaft in the SFB/TR 09. 

A Finite volume simulations 

For the matching in a finite volume, we performed one set of simulations of (quenched) 
QCD and one of HQET. In the case of the relativistic theory, we used L = L\, defined 
by </ 2 (Li/4) = 1.8811 The parameters of these simulations have been taken from [31] 
(see Table [3]). The difference is that here L = L\ = 2Lq (and T = L/2 and T = L/2±a 
in addition to T = L) compared to L = Lq in [31]. 

9 Zi differs slightly from L\ defined in the main text by L\ — 0.718ro. This mismatch is however 
corrected, as explained later in this appendix and in Add. [Pi 
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0.124528 
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7.6547 


0.133632 


1.8811(28) 


0.6713(8) 


-0.622 


1.0916 


0.127991 
0.126967 
0.126222 


40 


7.8439 


0.133373 


1.8811(22) 


0.6679(8) 


-0.623 


1.0900 


0.128989 
0.128214 
0.127656 



Table 3: Bare parameters used in the computation of the QCD observables for L = L\ 



The parameters for the resolution L\ja = 20 cannot be found in the mentioned 
reference. For that point, the gauge coupling f5 has been chosen such that 5 2 (Li/4) = 
1.8811 for L\/Aa = 5, see [20]. The renormalization constant Zp and n\ = k c have 
been computed here, while b m and Z have been extrapolated from the values in Table 2 
of [31]. These factors are put into the relationship between the bare mass m q h and the 
RGI mass [20,63], 

M = hZ m m qih (1 + 6 m am q>h ) , (A.l) 



where 



z z A 
z P 



and am. 



q,h 



1 

2 V k c 



(A.2) 



The renormalization constant Z\(qq) is known non-perturbatively from [30], while 

M 



h 



m(/i ) 



1.544(14) , no = 2/Li 



(A.3) 



relates the running quark mass in the Schrodinger functional scheme [20] at the scale 
^0) t° the renormalization group invariant quark mass M0 

For all values of L\/a three hopping parameters Kh have then been fixed in order 
to achieve 

z = L l M = 10.4, 12.1, 13.3. (A.4) 

We collect these parameters in Table [3l whereas the results for the quantities needed 
in the matching step are summarized in Tables H] and [5j The errors there include system- 
atics due to the uncertainties in the Z-factors, in particular, the error on the universal 
factor h has been propagated only after performing the continuum limit extrapolations. 



10 In h = M/m(fj,o) we take the small difference between the above defined Lo and the value Lo = Li/2 
into account. It causes a change of less than 1% of the value of h used in [31]. 
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Table 4: Simulation results of the finite volume (L = L\) relativistic observables needed 
in our main strategy. The continuum limits, obtained by linear extrapolation in (a/L) 2 
of the results for L/a > 24, are indicated by CL. 



L/a 


z 










$ 3 








9i = 


01 = 1/2 


0! = 1 


0i = 


0i = 1/2 


0i = 1 






02 = 1/2 


02 = 1 


02 = 








20 


10.4 


0.1699(9) 


0.4299(12) 


-0.5998(20) 


8.059(37) 


8.293(37) 


8.993(37) 


20 


12.1 


0.1668(9) 


0.4198(11) 


-0.5867(20) 


9.315(37) 


9.545(37) 


10.234(37) 


20 


13.3 


0.1649(9) 


0.4137(11) 


-0.5787(19) 


10.271(37) 


10.500(37) 


11.180(37) 


24 


10.4 


0.1739(23) 


0.4391(31) 


-0.6130(54) 


7.864(39) 


8.102(38) 


8.822(39) 


24 


12.1 


0.1710(23) 


0.4295(30) 


-0.6005(52) 


9.027(39) 


9.263(38) 


9.971(39) 


24 


13.3 


0.1693(22) 


0.4239(29) 


-0.5931(51) 


9.874(39) 


10.109(38) 


10.809(38) 


32 


10.4 


0.1760(41) 


0.4494(48) 


-0.6254(90) 


7.713(43) 


7.941(41) 


8.661(42) 


32 


12.1 


0.1733(40) 


0.4403(46) 


-0.6135(87) 


8.818(42) 


9.045(41) 


9.753(42) 


32 


13.3 


0.1717(40) 


0.4349(45) 


-0.6066(85) 


9.603(42) 


9.828(41) 


10.531(42) 


40 


10.4 


0.1790(70) 


0.4493(72) 


-0.6283(142) 


7.656(45) 


7.894(42) 


8.624(44) 


40 


12.1 


0.1763(68) 


0.4403(70) 


-0.6166(138) 


8.743(45) 


8.979(42) 


9.698(44) 


40 


13.3 


0.1747(67) 


0.4352(68) 


-0.6099(136) 


9.509(45) 


9.744(42) 


10.456(44) 


CL 


10.4 


0.1801(75) 


0.4587(84) 


-0.6392(159) 


7.533(89) 


7.765(86) 


8.496(88) 


CL 


12.1 


0.1776(73) 


0.4502(81) 


-0.6280(154) 


8.573(91) 


8.805(88) 


9.524(89) 


CL 


13.3 


0.1761(72) 


0.4452(79) 


-0.6218(151) 


9.289(93) 


9.519(91) 


10.234(92) 



Table 5: Same as Table U] in the case of the alternative strategy. 
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L/a 








^»stat 








0i 


= 


0i = 


= 1/2 


0i 


= 1 




02 = 


= 1/2 


02 


= 1 


02 


= 




HYP1 


HYP2 


HYP1 


HYP2 


HYP1 


HYP2 


6 


0.06936(5) 


0.06939(4) 


0.18583(7) 


0.18591(7) 


-0.25519(12) 


-0.25530(11) 


8 


0.07572(6) 


0.07574(6) 


0.20452(11) 


0.20457(11) 


-0.28024(17) 


-0.28031(17) 


10 


0.07821(5) 


0.07822(5) 


0.21246(8) 


0.21249(8) 


-0.29067(13) 


-0.29071(13) 


12 


0.07934(8) 


0.07935(8) 


0.21622(13) 


0.21625(13) 


-0.29556(21) 


-0.29559(20) 


CL 


0.08238(12) 


0.22596(21) 


-0.30835(32) 



Table 6: Lattice results of Rf 1 for L = L\. The continuum limits are obtained by a 
linear extrapolation in (a/L) 2 of the results for L/a > 8. 



LI a 






pstat 








0i 


= 


0i = 1/2 


0i 


===== 1 




02 = 


= 1/2 


02 = 1 


02 


= 




HYP1 


HYP2 


HYP1 HYP 2 


HYP1 


HYP 2 


6 


0.1502(3) 


0.1543(3) 


0.3562(3) 0.3688(3) 


-0.5231(6) 


-0.5231(6) 


8 


0.1544(4) 


0.1575(4) 


0.3672(4) 0.3765(4) 


-0.5216(7) 


-0.5340(8) 


10 


0.1571(5) 


0.1595(5) 


0.3724(6) 0.3710(6) 


-0.5295(10) 


-0.5391(10) 


12 


0.1561(8) 


0.1579(8) 


0.3729(8) 0.3786(9) 


-0.5289(15) 


-0.5365(16) 


CL 


0.1606(6) 


0.3827(6) 


-0.5432(11) 



Table 7: Lattice results of i?~ . The details are the same as in Table El 



Ensembles of roughly 2000 (for L\/a = 20) to few hundreds (for L\/a = 40) gauge con- 
figurations have been generated for this part of the computation. The lattice L\/a = 20 
is not used in the extrapolations but rather to check for the smallness of higher order 
cutoff effects for L\/a > 24. 

Concerning the simulation of HQET, we have computed the various quantities in 
the two required volumes. The first one, where we match the effective theory with 
QCD, has a space extent L\. The second one is such that L2 = 2L\. The value of the 
Schrodinger functional renormalized coupling is fixed at g 2 {Li) = 3.48, and we have 
used the resolutions L\ja = 6, 8, 10, 12. The corresponding values of (3 as well 
can be found in Table A.l of [1]. All these quantities are computed with two different 
actions, HYP1 and HYP2. The continuum values are then obtained by constraining the 
fits to give the same values for these actions. We note that the results for HYP1 and 
HYP2 are statistically correlated. 

For the computation of the step scaling functions one uses the same (5, n and Lija = 
2L\/a. All these computations are done with several thousand gauge configurations. 
Note that, even if L\ is the same in QCD and in HQET, the typical lattice spacings 
are much larger in the effective theory. The results of i?| tat and i? av can be found in 
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L/a 


E m (3.48, a/L) 


HYP1 HYP 2 


8 


0.431(11) 0.411(11) 


10 


0.437(11) 0.424(10) 


12 


0.422(16) 0.418(16) 


CL 


0.430(25) 



Table 8: Lattice results of the step scaling function E m . The bare parameters are 
described in the text. The continuum limit is obtained by a linear extrapolation in 
(a/L) 2 of the results for L/a > 8. 



L/a 






£f n (3.48, a/L) 








01 


= 


0i = 1/2 


0i 


= 1 




02 = 


= 1/2 


02 = 1 


02 


= 




HYP1 


HYP2 


HYP1 HYP 2 


HYP1 


HYP 2 


6 


0.6241(17) 


0.6245(11) 


0.6219(60) 0.6223(5) 


0.6225(8) 


0.6228(6) 


8 


0.5790(20) 


0.5797(13) 


0.5789(65) 0.5793(5) 


0.5789(10) 


0.5794(7) 


10 


0.5587(47) 


0.5586(22) 


0.5585(14) 0.5588(9) 


0.5586(22) 


0.5590(14) 


12 


0.5364(66) 


0.5342(39) 


0.5424(19) 0.5417(12) 


0.5409(30) 


0.5398(18) 


CL 


0.457(10) 


0.471(3) 


0.467(5) 



Table 9: Lattice results of the step scaling function S^ in . The continuum limits are 
obtained by a linear extrapolation in a/L of the results for L/a > 8. 



L/a 




£^(3.48, a/L) 








0i = 


0i = 1/2 


0i 


= 1 




02 = 1/2 


02 = 1 


02 


= 




HYP1 HYP2 


HYP1 HYP2 


HYP1 


HYP2 


8 


4.81(44) 4.72(32) 


1.58(15) 1.55(10) 


-1.19(11) 


-1.17(8) 


10 


4.34(58) 4.20(39) 


1.43(19) 1.39(13) 


-1.08(15) 


-1.04(10) 


12 


4.79(86) 3.98(58) 


1.58(28) 1.31(19) 


-1.19(21) 


-0.99(14) 


CL 


2.9(1.5) 


0.96(50) 


-0.71(38) 



Table 10: Same as Table M for 



Tables [6] and [7J The values of the step scaling functions are collected in Tables [HJ [9] 
andlOl 

Finally there are simulations in small volume to obtain the subtractions r® (L2) 
and T\ m (L2). These are done with L2 = 1.436 vq and j3 determined from the knowledge 
of r^/a [48]. The parameters, including listed in Table 6 of [27]. The values 

of (3 do of course agree with the ones employed in the large volume, which we describe 
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13 




L 3 


x T 




a 2 


£"kin 










HYP1 


HYP 2 


HYP1 


HYP 2 


6.0219 


0.133849 


16 3 x 


24 [32] 


0.4345(21) 


0.4029(32) 


0.750(4) 


0.774(3) 


6.0219 


0.133849 


24 3 


x 32 


0.4378(25) 


0.4034(20) 


0.746(7) 


0.776(5) 


6.2885 


0.1349798 


24 3 


x 48 


0.3295(21) 


0.3034(29) 


0.643(7) 


0.676(5) 


6.4956 


0.1350299 


32 3 


x 64 


0.2724(20) 


0.2461(14) 


0.599(10) 


0.620(11) 



Table 11: Parameters of the large volume simulations. Where present, the numbers in 
brackets refer to a second dataset at the same (/?, k) values. 

in the next appendix. 

B Large volume simulations and extraction of matrix elements 

The parameters for the simulations in large volume are collected in Table [TT] together 
with the results for -E s tat and -Ekin- The lattice extension L/a and (3 are such that 
L = 4Li ~ 3/2 fm except for the second lattice where we have L = 6L\ ~ 2fm. This 
lattice is used only to check for the absence of finite size effects. We see from Table [TT] 
that finite size effects are indeed very small, the difference between the results from 
the L/a = 16 and the L/a = 24 lattices at f3 = 6.0219 is consistent with zero within 
at most one standard deviation (diktat from HYP1). The number of configurations 
generated ranges from 4300 at (3 = 6.0219 to 2200 at (3 = 6.4956 (for the larger volume 
at (3 = 6.0219 we had 1300 configurations). Since our phenomenological input is the 
mass of the (spin averaged) B s meson, we set k to k s in order to reproduce the quenched 
value of the strange quark mass from Ref. [64], i.e. 

M s r = 0.35(1) , (B.l) 

with M s the renormalization group invariant strange quark mass defined as in Ap- 
pendix [A] after replacing by k s . 

The numbers for -E s tat and £kin have been obtained by applying two different fitting 
procedures to two independent datasets (where available). The quoted errors are such 
that both the results are covered and they therefore provide a reasonable estimate of 
the systematics associated with the fits. We now sketch these procedures. 

Let us consider in QCD the effective "mass" T(xo) obtained from the correlation 
function F iw (xq) in eq. (jC.ip and its quantum-mechanical decomposition 

T(xo) = -^^Fav = E Q + Ae~ Axo + ... (B.2) 

where Eq is the energy of the ground state, A is the gap between the ground and the 
first excited states and the dots refer to contributions from higher states. The \/ra\ ) 
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expansion reads 

r(* ) = E Btat + a, kin £ kin + (A** + a, kin A kin ) e - AStata; ° (1 - ^ ki nX0 A kin ) + . . . 

= r stat (x )+^ kin r kin (xo) + ... (B.3) 

where r stat and r kin are denned in analogy to eqs. (|3. 101 [37TT]) in terms of the correlators 
/I tat (x ) and /l in (x ). 

In the correlation function fs a( x o) the same states contribute as in /a(xo)- Per- 
forming again first the quantum- mechanical decomposition and then the 1 /m^ expansion 
of these correlators, it is easy to see that the ratios 



psr[xo) = Ww) and ^ in(xo) = pr {xo) 

have the following form 



(B.4) 



P| tat = b x + 6 2 e- Astatxo , (B.5) 
pkm = 6 3 + 6 4 e- Astata;o - 6 2 A kin x e- Astatxo . (B.6) 

They can therefore be used to further constrain A stat and A km . We are thus lead to 
perform a combined fit 

T stat = 65 + bee-" 1 """ , (B.7) 
T kin = b 7 + 6 8 e- Astata:o - 6 6 A kin x e- Astata;o , (B.8) 

together with eq. (|B.5|) and (|B.6p . with non- linear parameters a± = A stat and 02 = A km 
and the linear parameters &«, which contain the desired 65 = E sta _ t and 67 = -E^in- 

Since the correction terms are nevertheless not so easy to compute at the smaller 
lattice spacings, we perform the above fit first at (3 = 6.0219 and extract aA stat and 
a 2 A kin . We then use that these quantities scale roughly (i.e. tq A stat ~constant and 
Tq A kin ~constant). To implement this, we input the scaled means as priors [65] in a 
second step where we add 

(prior \ 2 
cij — aX 

W - / (B-9) 

a. prior / j ,r prior \o ' \ / 

=1,2 (K 



tprior (X„P lior \2 



to the standard x 2 . The uncertainty 5a? nov is taken from the fit result at (3 = 6.0219. 
However, in order to remain on the safe side, it is not scaled but kept constant at the 
smaller lattice spacing. Thus 5a^ 10T /a^ lov oc I /a 2 for example. The constraint due to 
the priors becomes weaker as we approach the continuum. 

Here and in the following procedure the fit range is chosen to keep a minimum 
physical distance from the boundaries, namely x$ > i m i n « 2tq. The stability of the 
results is checked by varying i m i n to t m \ n — ro/2. As an example we show in figured] 
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the results for P| tat , Pg n , T stat and T km at (3 = 6.2885. One observes that Pf at , T stat 
provide very good constraints of the parameters A stat , 62 ; b%. The remaining ones are 
then effectively linear fit parameters. Nevertheless, the error band of i^kin (dashed line) 
resulting from the fit is not that small. 

An alternative strategy is used to get a second estimate of EUat at the two coarser 
lattice spacings, where we have two independent datasets. Exploiting again the remark 
before eq. (jB.4[) we construct an effective mass T$a from the correlator f$ a(^o) in the 
very same way as r stat is obtained from /i (#o)- The idea is to combine the two 
effective masses in order to eliminate the contribution from the first excited state and 
then perform a fit to a constant (in the mentioned fit range). In practice we minimize 
the quantity 

^max 

Q = [r stat (*o) + aT 5A (x ) - B] 2 , (B.10) 

^•0 ^min 

with respect to a, B. Finally the weighted average of [r stat (xo) + oT^a^o)] /(l + a) 
yields the estimate of E staX . The quality of the result is comparable to that obtained in 
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the first approach. 



C Alternative strategy 

We briefly introduce an alternative strategy, based on the correlation functions /a, ky 
in addition to fi,k±. With 

F av (x , 9) = - A log (- [/ A ] R (so, 0) x [£v]r (xo, 0)) (C.l) 

we introduce 

R av (L,e 1 ,9 2 ) = F av (x ,#i) -F av (x ,# 2 ) at x = L/2, T = L (C.2) 

r av (L,# ) = -^^F av (x ,#o) at x = L/2, T = L. (C.3) 

Keeping 3>i, from the standard strategy, we define the set of observables 

5i(L,M b ) = $i(L,M b ), (C.4) 

$ 2 (L,M b ) = R av (L,e 1 ,e 2 )-Rt at (L,e 1 ,e 2 ), (c.5) 

<D 3 (L,M b ) = Lr av (L,0 o ), (C.6) 
with the l/m b expansion 

$ 2 (L,M b ) = Wkin i? A in (L^ 1 ,0 2 ) + Ca i Q ET i? 5A (i^i^2) (C.7) 
$ 3 (L, M b ) = L [m barc + T stat (L, ) + w kin r kin (L, O ) + c a i v QET r M (L, 9 )] ,(C8) 

where due to the spin average the combination 

c^ ET = \lc R A QET + 3c^] (C.9) 

is present. The so far undefined terms R^, R^ 11 , r km , Rsa, T<5A are straightforwardly 
obtained from our definitions. 

The alternative observables change from L to 2L via 

$,(2L,M b ) = Y,aij{u)*j{L,M h ) + 5= m («) . ( C - 10 ) 

(Jij(u) = lim T,ij(u,a/L) (Cll) 

with the step scaling functions (we drop arguments 6>i, 6> 2 and u = g 2 (L) is understood) 
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do 


r M<°> 










roM^ 








0i = 


01 = 1/2 


X = 1 


0i = 


0i = 1/2 


0! =1 






e 2 = 1/2 


02 = 1 


02 =0 


2 = 1/2 


02 = 1 


02 =0 





17.05(25) 


0.17(6) 


0.17(6) 


0.17(6) 


0.02(9) 


0.02(8) 


0.02(9) 


1/2 


17.01(22) 


0.20(7) 


0.18(6) 


0.19(7) 


0.02(10) 


0.02(9) 


0.02(9) 


1 


16.78(28) 


0.34(11) 


0.30(7) 


0.32(8) 


0.06(12) 


0.06(9) 


0.06(10) 



Table 12: RGI results of Mb i n the static approximation and of the l/m b correction 
for the alternative strategy. 



Su(u,o/I) 

E 2 i(w,o/-L) 
£ 22 (u,a/L) 

S 3 i(n,a/L) 

£ 32 (u,a/L) 
£ 33 (u,a/L) 

<7 m («) 



Rf n {2L)/Rf n {L) = Sf n (u,a/L) 



jkin 



kin/ 



1 



{i?| in (2L)-i?l in (L) Z 22 (u,a/L)} 



i? kin (L) 
R SA (2L)/R SA (L) 

2L{r kin (2L) - T kin (L)} 
fl kin (L) 

r 5A (2L)-r 5A (L) 



s 32 ( n , a / L )<!^l 

V ' J R\ m (L) 



2L 
2 



lim 2L [r stat (2L) - r stat (L)l . 

a/L^0 



The final relation for the B-meson mass is eq. (|2.30p with 



L 2 m^\M h ) 
L 2 mg h) (M h ) 
L 2 m^ a) {M h ) 

L 2 mg 6) (M b ) 



a m (u 1 ) + 2<S> 3 (L 1 ,M h ), 
L 2 [E stat -T stat (L 2 )], 

(T3l(ui)*l(Ll,M b ) + f73 2 (ui)$2(il,Mb) 



(C.12) 

(C.13) 
(C.14) 

(C.15) 

(C.16) 
(C.17) 
(C.18) 



(C.19) 
(C.20) 
(C.21) 



L 2 



E km - r kin (L 2 _) r,5 A (L 2 ) R A m (L 2 ) 



-L. 



Rf n (L 2 
Tsa(L 2 
Rsa(L 2 



a kin (m)$i(Li,M b ) 



R SA (L 2 )R^(L 2 )_ 
C j 2 iK)i 1 (L 1 ,M b )+ C T 22 (n 1 )$ 2 (L 1 ,M b ) . (C.22) 



Although the results have been already given in Table O the reader will find more 
details in Table 12 

D Propagating uncertainties in Li/r and g 2 [Li) 

In our simulations we have fixed L\ by g 2 {L\/4) = 1.8811, because the corresponding 
bare parameters j3, k are available in the literature. We here give the estimate of the 
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small effect caused by L\ ^ L\ in the static approximation. From the polynomial 
interpolations of the step scaling function of the coupling, cr(u) [20], we estimate the 
corresponding mismatch in couplings as 

u-u = 5 2 (Li) - g 2 {L{) = ct(<t(1.8811)) - 3.48 = -0.17(5) . (D.l) 

Let us write 

— = p(u,z)[l+K(u)} at u = u (D.2) 
m B 



with 



zri \ r| tat (Li) — -Estat i \ Z IT\^\ 

K{u) = , p(u,z) = — -. (D.3) 



The relation -r- — = gives 

cm me ° 

p{u, z) du wb dn 

Denoting by AMb the correction we have to add to when it is computed with u ^ u 
(as we did), we get from the above equations 

— AM b = [u-u] X p(u)K'(u) , (D.5) 
m B 1 

where K'(u) is easily estimated by taking a numerical derivative of T| tat . From the 
difference of L/a = 12 and L/a = 10 at fixed g$ (with g 2 \L/ a =i2 = 3.48) and with 
p(u, z) ~ 1.44 we arrive at the small shift 

r AM b = -0.055(17). (D.6) 

A similar error is be taken into account due to the 2% uncertainty in the relation 
L,2 = 1.436ro [48]. In the same way it leads to a statistical error of 

r AM b = 0.016. (D.7) 

The two contributions eq. (|D.6p . eq. (|D.7p are combined to 

r AM b = -0.055(23), (D.8) 



which we have taken into account in Sect. 15.31 Because of the smallnes of these effects, 
they can be neglected in the 1/mb-corrections. 

In the case of our alternative strategy, the shift depdends on the value of #o- We 

find 

O = r AM b = -0.042(20) , (D.9) 

9 = 1/2 r AM b = 0.009(11) , (D.10) 

6 = 1 r AM h = 0.150(45) . (D.ll) 
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